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1. INTRODUCTION 

Electric power, worldwide, has become a consumer good for both everyday life and the economy of 
the countries, and the smallest problems of electrical origin have a great influence on the continuity of 
economic activities. Therefore, the possession of reliable and economical electricity networks operating 
properly and ensuring a continuity of service and a good quality of energy has become essential in order to 
contribute positively to the development of our modern societies. 

The distribution network must evolve towards a flexible and intelligent network that best integrates 
local and / or renewable energies. The opening up of the electricity market and the growing environmental 
concerns linked to global climate change are leading to significant changes, particularly in distribution 
networks, with the massive influx of decentralized products. This evolution can be envisaged by developing 
intelligent systems, capable of minimizing the impacts generated by the insertion of decentralized 
productions and / or by the search for new architectures. These two solutions should allow the increase of the 
decentralized production rate in the distribution network under the best economic and security conditions. 

The distribution network was not originally designed to house production units, but to provide 
unidirectional electricity from the distribution network to medium and low voltage consumers. The use of 
Distributed intelligent systems alone will not be enough to eliminate all the problems that the distribution 
network will face in case of significant penetration of decentralized productions [1]. 

Distribution networks are branched radial topology networks and have several load nodes. The problem that 
arises is how to choose the best location of GED (decentralized energy generator) in a distribution network. 
This problem is the subject of research thematic of this article, for this we have used the software for 
calculating the power flow and calculate some parameters to determine the best node of GED connection. 
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Several methods have been used to choose the optimal location of the wind turbine in the distribution 
network in order to reduce the active losses. [2], Try to provide a comprehensive approach to solve the 
problem of placement and determining the capacity of wind units in the network. [3], presents a voltage 
sensitivity analysis with respect to the real power injected with renewable energies to determine the optimal 
integration of distributed generation (DG) in distribution systems (DS). 

In this article, a technique used is based on the particle swarm optimization method (PSO) to 
determine an optimal allocation of wind-based distributed generators in order to reduce the real power loss of 
the distribution system. The developed method is implemented on the IEEE 33 bus test system. The work 
plan for this article has been organized as follows: The basic mathematical model of radial distribution 
systems is presented in the Section 2. Wind generation system modeling is given in Section 3. In section 4 
Loss sensitivity factor and PSO method is used for optimal allocation of wind. Results and Discussion 
explained in sectionS. Finally, concludes this paper and proposes perspectives for continuation work. 


2. Mathematical model 
2.1. Power flow formulation 

The goal of this article is to find the optimal location of the wind turbine in the network in order to 
minimize active losses. The technique proposed for calculating the power flow in the radial distribution 
network is based mainly on knowledge of the topology or architecture of the network. The single-line 
diagram is shown in Figure 1. The power flow equations are given by [4-6]. 
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Figure 1. Single line diagram of radial feeder 
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Where , P; and Q; are the real and reactive power flowing out of bus i; P,;4, and Q,;4, are the real 
and reactive load powers at bus i + 1.line section between buses i and i + 1has resistance R; ;,,and 
reactance X; ;4,.voltage magnitude of bus i is |V;|.for convergence of power flow ,the power balance (1)and 
(2) must be satisfied in addition, magnitudes of sending and receiving end bus voltages must satisfy (3). 

The well known generalized formulas for real and reactive power loss in the line section between buses i and 
i + 1 are calculated by using the following equations. 
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The total active and reactive losses in the distribution network can be easily found by adding all the 
power losses in all branches and they are expressed as follows. 


Prtoss = a Pross (i,i T 1) (6) 
Qtioss = ae Qross (i, i+ 1) (7) 
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Where, Pios5 1s the total real power loss and Qioss 1s the loss incurred due to reactive power throughout the 
network. In this study we consider wind supplying real power with unity power factor. 


2.2. Electrical grid operating constraints 


To ensure the stability of the network and the safety of equipment, voltage must be maintained within 
limits set by the rules of operation of the power grids. 


Vmin < Vi < Vmax (8) 


Where Vmin and Vmax are the minimum and maximum values of bus voltage amplitudes, respectively. 
And the power transited in a line must not, under any circumstances, exceed the maximum limit. 


Ss (9) 
With 
— (p2 2x4 1/ 
Sij = (Pg + QG) /2 (10) 
Where 
S;; : Apparent power transmitted in line iż - j. 
And 


Si; : Maximum apparent power transited in line i - j. 


3. WIND GENERATION SYSTEM MODELLING 

The power generated by the wind turbine depends on their type and location, as well as the wind 
speed [7-12]. The electrical output power of a wind turbine is determined by using (11). 
With Veins Veour Vy are different speeds corresponding to the wind turbine and P,.,;.qis rated out pout power 
of wind turbine and it can be determined from (12). 


OV 5 < Vee, OF Vig > Vout 

Vw-Vcin 

Pw = err, ee Voi < Vy < Vy (11) 
Prated Vy < W 2 Veout 





Pratea = 0.5pAV3C, 02) 


4. LOSS SENSITIVITY FACTOR 
The Active power losses in electrical networks are given by (12). 


P, = YL Dai (PP; + Q:Q;) + Bij (Qi P; + P:Q;)| (13) 


Where, 


r 


ij — sS. 
yv; °S (Ôi ô) (14) 





Qij = 


Bij = ysin (5;-6;) (15) 


And 
Zij = Tij + Xij Are the ijt” of [Zous| matrix 
P; = Pej — Ppi And Q; = Qei — Qoi 

Pg & Qg; : Active and reactive power of the generator. 


Ppi & Qp; : Active and reactive load power 
The sensitivity factor of real power loss with respect to real power injection from the DG is given by 
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ap 
aj = T = 2; P; + 2 Dj=1(aijP; — BiQ) (16) 


j+i 


The sensitivity factor is calculated for each node using the values obtained from the calculation of the basic 
flow. The node with the lowest loss sensitivity factor will be the best node for DG placement [13-15]. 


5. PARTICLE SWARM OPTIMIZATION 

5.1. Introduction 

The PSO algorithm is initialized by a population of random potential solutions, interpreted as 
particles moving in the research space. Each particle is attracted to its best position discovered in the past as 
well as to the best position discovered by the particles of its vicinity (or of the whole swarm, in the global 
version of the algorithm). The PSO algorithm includes several adjustment parameters which act on the 
exploration-exploitation compromise [16-22]. Exploration is the ability to test different regions of space in 
search of good candidate solutions. Exploitation is the ability to focus research around promising solutions in 
order to get as close to the optimum as possible. The choice of parameters remains largely empirical. A 
complete analysis of the algorithm was made by Kennedy 
In general, a PSO algorithm has three main steps: 

e First step: consists in authorizing all the particles arbitrarily distributed across the search space. 

e Second step: consists in evaluating the fitness value of the particles in order to determine the best 
position of each particle and to reveal the particle which has the best overall fitness value in the current 
swarm. 

e Third step: updates the speed of the particles which are then used as information to change the position 
of all the particles. 

In PSO, two different definitions are used: the individual best and the global best. As a particle 
moves through the search space, it compares its fitness value at the current position to the best fitness value it 
has ever attained previously. Adjustment of the position and speed of each particle is calculated using the 
current positionx,;, the best particle position so far Pbest and the best overall particle position in the 
populationGbest. Velocity of each particle in the next generation can be calculated as: 


ie = wV‘ + cn (Pbest — x )+ cr, (Ghest — X ) 


(17) 


With 7, and r, are two numbers between 0 and | while c} and c, are positive constants, respectively. 


Inertia weight w is is formulated by the following expression: 


W — Wmi 
W = Wax — oe iter 
Mer nax (l 8) 
Where w nax and w,,:, are maximum and minimum of inertia weight, iter nax and iter are maximum iteration 


number and current iteration, respectively. 


5.2. PSO Procedure 

The calculation procedure for the PSO calculation algorithm is as follows. 
Step 1: Introduce data for all nodes and branches as well as the limits of nodal voltages. 
Step 2: Calculating the power flow in the distribution network and determining total active losses. 
Step 3: generate an initial population of particles with random positions and velocity. Set the iteration index 
t= 0. 
Step 4: For each particle, update the bus data (for wind turbine) based on its locations and setting values, if 
all nodal voltages are within the permissible limits, estimate the total loss. 
Step 5: Compare the fitness value of each particle with the personal best, Phest . If the fitness value is lower 
than Pbest, set this value as the current Pbest, and record the particle position corresponding to this Pbest 
value. 
Step 6: Select the minimum value of Pbest from all particles to be the current global best, Gbest , and record 
the particle position corresponding to this Gbest value. 
Step 7: Update the velocity and position of particle. 


Optimal placement of wind turbine in a radial distribution network using PSO method (Yahiaoui Merzoug) 


1078 O ISSN: 2088-8694 


Step 8: If the maximum number of iterations is reached, the particle associated with the current Gbest is the 
optimal solution and then go to Step 9. Otherwise, set t=t+/ and return to Step 4. 
Step 9: Displayed the best position for each particle. 


STAR 
Input system data (1) 


Calculate the loss using Distribution load flow (2) 





Initialise particule population (3) 


Calculate the total loss (4) : 
Record Pbest (5), Gbest(6) : 
Update particle position and velocity (7) : 


No : 
Stop test (8) = 


. ei ee Ue Yes 
Particle Swarm ptimization 
Print out location of wind (9) 


Figure 2: PSO-Optimization Wind computational procedure 





6. RESULTS AND SIMULATIONS 
Optimization of the location of the wind is done by the PSO algorithm, always programmed under 


MATLAB and tested on the IEEE33-bus distribution network system containing 33 buses and 32 branches as 
shown in Figure 3. The bus and line data are taken from [23-24]. The real and reactive power loads and 
voltage of radial distribution system are 3.72MW, 2.3Mvar and 12.66KV respectively. This program 
searches for the best wind turbine location on the network to minimize total active losses. 
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Figure 3. Single line diagram of 33 Bus electric distribution systems. 


Voltage sensitivity index was calculated at all nodes. Bus 18 was found to have the least VSI. Hence 
Wind Turbine was placed at this bus. The active power losses before the placement (basic calculation) of the 
distributed generator (wind) is 208,492 KW and the minimum voltage is 0.91075 pu at node level number 18. 
The first step in determining the wind power injected into the grid using the (11), with an average annual 
wind speed considered to be 6.02 m/s to determine the optimal location. [25]. 
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Figure 5. Bus Voltages Before and After Optimal 


Figure 4. Loss Sensitivity at all nodes DG Installation in the 33-Bus System 


After insertion of the wind, it is noted that there is an improvement in all the voltages and an active loss 
reduction of the 33.355% as shown in the following table 1. 


Table 1. Power loss reduction. 
Without wind With wind 
Power loss (KW) 208.4592 138.9275 
Power loss reduction - 33.355 % 


Minimum voltage 0.91075 


The comparison for the real power losses in the branches is shown in Figure 6 for the initial network 
(with out wind) and with wind (case 1 and case 2), respectively, The power losses in almost every branch in 
case 2 reduced, except at 18, 19, 20, 21, 33, 34 and 35, where there was a small increase in losses. 
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Figure 6. Real power losses in 33-bus system 


6. CONCLUSION 

The optimal insertion of wind turbines into the distribution network, the technique used to 
contribute to the management of distribution networks with the aim of minimizing total power losses by 
ensuring a good voltage profile, are carried out in the process of based on optimization by a meta heuristic 
method called pso algorithm. Many perspectives open up on the subject of the optimal insertion of 
decentralized production in the distribution network, in particular with regard to the types of DGs installed 
(photovoltaic, wind, etc.) and their maximum numbers to be installed, thus reconfiguring the distribution 
network taking into account other constraints such as transit limits and the assurance of continuity of service 
in the event of faults. 
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